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Abstract- Planar laser-induced fluorescence (PLIF) is 
demonstrated to be a valuable tool for studying the onset of 
transition to turbulence. For this study, we have used PLIF of 
nitric oxide (NO) to image under expanded axisymmetric free jets 
issuing into a low-pressure chamber through a smooth 
converging nozzle with a sonic orifice. Flows were studied over a 
range of Reynolds numbers and nozzle-exit-to-ambient pressure 
ratios with the aim of empirically determining criteria governing 
the onset of turbulence. We have developed an image processing 
technique, involving calculation of the standard deviation of the 
intensity in PLIF images, in order to aid in the identification of 
turbulence. We have used the resulting images to identify 
laminar, transitional and turbulent flow regimes. Jet scaling 
parameters were used to define a rescaled Reynolds number that 
incorporates the influence of a varying pressure ratio. An 
empirical correlation was found between transition length and 
this rescaled Reynolds number for highly underexpanded jets. 

I. Introduction 

Highly underexpanded axisymmetric free jets are important 
flows of interest for fundamental fluid mechanical studies as 
well as for practical applications. Of particular importance is 
an understanding of the development of instabilities in these 
flows. While the literature contains extensive studies of the 
transition to turbulence in boundary layers of flat plate flows, 
the literature containing experimental studies of transition in 
highly underexpanded axisymmetric free jets is considerably 
more sparse [1-6]. Understanding the stability characteristics 
of such flows became the motivation for this work, conducted 
as part of the NASA Shuttle Return to Flight (RTF) 
investigation into the effects of flow through vehicle wing 
breaches. 

While CFD (computational fluid dynamics) codes can 
simulate turbulent flow, accurately predicting the transition 
from laminar to turbulent flow remains an unsolved problem. 
The empirical correlation resulting from these tests can 
therefore potentially be used to determine where turbulence 
models should be implemented in CFD simulations as well as 
providing data for comparison with CFD results. 

Flow visualization was required for these tests. The flow 
environments encountered in these tests include regions of low 
static pressure, turbulent and/or three-dimensional flow 
structures, and regions of interest that lack strong density 
gradients. Such conditions, though frequently encountered in 


aerospace simulation facilities, cannot be satisfactorily 
visualized using traditional path-averaged techniques such as 
schlieren and shadowgraph, which rely on sufficiently strong 
density gradients. An alternative approach is therefore 
required to visualize these types of flows. Planar laser-induced 
fluorescence (PLIF) is a powerful flow visualization technique 
that provides a means of making non-intrusive measurements 
with sub-millimeter spatial resolution and flow-stopping 
temporal resolution in many of these challenging testing 
regimes. [7] 


II. PLIF Measurement System 


A. Test Apparatus 

Tests were conducted at NASA Langley Research Center 
using the test section of the 15 -Inch Mach 6 Wind Tunnel as a 
vacuum chamber. Gas bottles were stored in a ventilated, 
alarmed toxic gas cabinet. Nitrogen (N 2 ) seeded with 0.5% 
nitric oxide (NO) was plumbed into a heated @500 K) 
stainless steel plenum, through a converging nozzle with a 
nominal exit Mach number of 1, into the vacuum chamber. 
Mass flow controllers adjusted the flow rates. A schematic 
including these elements of the test apparatus is shown in Fig. 
1 . 

Two parameters were varied: Reynolds number ( Re ex u ) and 
jet pressure ratio (. JPR ). Re ex it was defined in terms of the 
nozzle exit diameter, D e , as well as the density p e , velocity V e , 
and dynamic viscosity p e at the nozzle exit, as given by (1): 
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(Note that in keeping with conventions found in the literature, 
data are plotted in terms of a unit Reynolds number, with units 
of m" 1 . The unit exit Reynolds number is simply Re exit / D e . 
The nozzle used in these tests had an exit diameter of D e = 2.4 
mm.) Re ex it was varied by changing the mass flow rates. 

JPR was defined as the ratio of the static pressure at the 
nozzle exit, p e , to the ambient pressure in the vacuum 
chamber,^, according to (2): 

JPR = Ei- . (2) 
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Figure 1. Schematic of experimental setup. Gas is plumbed through a 
heated plenum and nozzle into a vacuum chamber. A laser sheet enters the 
top of the vacuum chamber and excites nitric oxide molecules in the flow. 
An intensified CCD camera positioned at right angles to the laser sheet 
images the fluorescence. 


JPR was varied by changing the pressure in the vacuum 
chamber (p a ). For a given Re exit , the maximum attainable JPR 
was limited by the minimum attainable pressure in the vacuum 
chamber, which was about 130 Pa (1 torr). 

A thermocouple and a pressure transducer measured the gas 
temperature T 0 and pressure p 0 in the plenum, upstream of the 
converging nozzle. Measuring these quantities made it 
possible to calculate conditions at the nozzle exit, and 
therefore the exit Reynolds number. Nozzle exit temperature 
T e (K) and velocity V e (m/s) can be calculated from these 
measured quantities by assuming inviscid adiabatic flow and 
an exit Mach number of one, according to (3) and (4) [8]: 
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V e =p^RT e , (4) 


where y is the ratio of specific heats of the gas (y = 1.4 for 
nitrogen), and R is the specific gas constant with a value of 
287 J/(kg-K). Furthermore, with the added assumption of 
isentropic flow, nozzle exit pressure p e (Pa) and density p e 
(kg/m 3 ) can be calculated using (5) and (6) [8]: 
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The dynamic viscosity ju e (Pa-s) at the nozzle exit was 
calculated using the power law approximation given in (7): 
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This power law accounts for the effect of temperature on 
dynamic viscosity while neglecting the slight dependence on 
pressure. Reference [9] gives an empirical value of n = 0.67 
for nitrogen, where T re f= 273 K and ju re f = 1.663 x 10" 5 Pa-s. 
Values for nitrogen were used because the gas consisted of 
99.5% N 2 , and because the NO that made up the remaining 
0.5% has a viscosity similar to that of N 2 [9]. 

B. Laser System 

The laser system consists of three main components: a 
pump laser, a pulsed dye laser (PDL), and a wavelength 
extender (WEX). An injection-seeded Nd:YAG laser 
operating at 10 Hz pumps a tunable PDL. The output of the 
dye laser and the residual infrared from the Nd:YAG are 
combined in a WEX containing both a doubling and a mixing 
crystal. The resulting output is tuned to a wavelength of 
226.256 nm, chosen to excite the strongly fluorescing spectral 
lines of NO near the Qi branch head. In this notation, Q 
indicates that the rotational quantum number J is the same for 
both the upper and lower states of the transition; i.e., A J=0. 
The subscript 1 indicates that, in both the upper and lower 
states involved in the transition, J is related to Nby J=N+l/2 , 
where N is the quantum number for total angular momentum, 
excluding spin. 

A monitoring gas cell system is used to ensure that the laser 
is tuned to the correct spectral line of NO. The gas cell 
contains a low-pressure mixture of 5% NO in N 2 . A quartz 
window serves as a beam splitter and sends a small portion of 
the laser energy through windows on either side of the gas 
cell. A photomultiplier tube (PMT) monitors the fluorescence 
intensity through a third window at right angles to the path of 
the laser beam. 

The components of this laser system are mounted within a 
two-level, enclosable, portable cart. A photograph of this 
portable PLIF system is shown in Fig. 2 with the panels 
removed to show the internal components. When all of the 
panels are in place, a single monochromatic ultraviolet laser 
beam exits the cart, creating a relatively safe operating 
environment. 

Optics form the beam into a 100 mm wide x ~1 mm thick 
laser sheet, which is directed vertically downward through a 
window in the top of the test section. Fluorescence is imaged 
onto a gated, intensified CCD at a viewing angle normal to the 
laser sheet. A 1 mm thick [Schott glass UG5] filter was 
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Figure 2. The portable PLIF system, shown with panels removed. 
Components include: (1) Nd:YAG laser; (2) dye circulators with Rhodamine 
590 and Rhodamine 610 laser dyes in a methanol solvent; (3) wavelength 
controller for the (4) pulsed dye laser; (5) wavelength extender; and (6) low- 
pressure monitoring gas cell. 

placed in front of the camera lens in order to attenuate 
scattered light at the laser’s frequency. The schematic of the 
experimental setup shown in Fig. 1 includes these elements as 
well. 

C. Image Acquisition 

Images were acquired at 10 Hz with a lqs camera gate and 
a spatial resolution of 7 pixels/mm. This temporal and spatial 
resolution was more than sufficient to resolve turbulent flow 
structures of interest. Further details of the laser and camera 
system can be found in [10]. 

Sets of 100 single-shot images were acquired for a range of 
unit Re ex it (0.07 x 10 6 m" 1 to 14.8 x 10 6 m" 1 ) and JPR (1.8 to 
38). Background images were also acquired on each day of 
testing for a range of vacuum chamber pressures. During the 
acquisition of these background images, the laser was firing 
but no gas was flowing through the nozzle. Any nonzero 
intensity in these background images is due to either camera 
dark current or the laser scatter and room light not blocked by 
the filter in front of the camera lens. Background images were 
created from the average of 100 single-shot images in order to 
smooth out random shot-to-shot variations in background 
intensity. 


D. PLIF Intensity 

The absolute intensity at any given pixel in a PLIF image 
has a functional dependence on many local flow parameters as 
given by (8): 

S f °c E p<$>Xm [ Z Ib B u Si ] , (8) 

where Sf is the fluorescence signal; E is the energy in a single 
laser pulse; p is the local gas density; O is the fluorescence 
yield; ^ N0 is the local mole fraction of NO; fg is the 
Boltzmann fraction; B n is the Einstein absorption coefficient 
for the probed transition; and gi is the spectral overlap integral. 
O, fg, and gi are functions of the local temperature; O and gi 
also depend on the local pressure. Since the fluorescence 
intensity depends on so many local flow parameters, variations 
in any one of these parameters will affect the resulting 
intensity in a PLIF image. While this means that the 
interpretation of PLIF images is not always straightforward, it 
also means that PLIF images are well suited to detecting 
fluctuations within a flowfield. 

E. Image Processing and Analysis 

Single-shot images were processed to correct for 
background scattered light and camera dark current as well as 
spatial variations in laser sheet intensity. Both the background 
image and the single-shot images were smoothed with a filter 
(a 3 pixels x 3 pixels rotationally symmetric low-pass 
Gaussian filter with a standard deviation of 1) prior to 
additional processing in order to reduce noise in the images. 
A background image was then subtracted from the filtered 
single-shot image. Images were subsequently corrected for 
laser nonuniformities. Conveniently, jet gas containing nitric 
oxide diffused relatively uniformly into the test chamber in 
regions away from the jet. The fluorescence from the diffuse 
nitric oxide in these regions provided a suitable laser-energy 
reference, allowing the spatial variation in the laser intensity 
to be corrected in the images. This was accomplished by 
selecting an area of the image above the core of the jet flow 
and finding the average intensity along each column in that 
region. (Note that raw images were 512 x 512 pixels; images 
presented in Figs. 3, 7, 8 and 9 have been cropped to show the 
regions of greatest interest.) Finally, the background- 
corrected images were divided by this laser sheet profile. 
Since the laser sheet is typically brightest in the center of the 
sheet, the result of dividing an image by the laser sheet profile 
was to increase the brightness of the left and right edges of the 
image relative to the central part of the image. 

After background subtraction and laser-sheet correction, 
each set of 100 images was processed to create an averaged 
image at that flow condition. Single-shot images were then 
compared to this averaged image. A map was created 
showing the standard deviation of the fluorescence intensity at 
each pixel in the images. This map was generated by first 
subtracting the averaged image from each single-shot image, 
squaring the difference at each pixel, averaging the squared 


differences of all 100 images, and finally, taking the square 
root of this average value at each pixel. The resulting map is 
then an indication of how much the intensity varied on 
average at a given location in the flow. Regions with small 
variations were interpreted to correspond with laminar flow, 
while the large variations in other regions were attributed to 
some degree of turbulence. 

III. Experimental Results 

Fig. 3 shows a typical image sequence. The top image is a 
single-shot image. Towards the right of the image, turbulent 
structures are clearly visible. The middle image in the 
sequence shows the average of 100 single-shot images taken at 
this same flow condition. The bottom image is the map 
showing the standard deviation of the intensity of 100 single- 
shot images at each pixel. 

For the purposes of this study, we chose to define laminar 
regions of the flow as those where fluctuations within the flow 
were comparable or less than the background fluctuations in 
regions outside the core flow of the jet. The onset of transition 
is then defined to occur at the downstream location where 
fluctuations begin to rise significantly above the background 
level. This is marked as x tram in Fig. 3. Finally, we identified 
the flow as being fully turbulent downstream of the location 
where instabilities have penetrated the core of the jet, as 
evidenced by an increase in the standard deviation at the jet 
centerline. This location is labeled x turb in Fig. 3. For jets that 
remained laminar within the field of view of the camera, the 
distance from the nozzle exit to the rightmost edge of the 
image was identified as xi am . 

For each test condition, the standard deviation image and 
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Figure 3. (top) Single-shot image of PLIF intensity; (middle) 100-shot 
average image for unit Re exit = 2.7 x 10 6 m' 1 and JPR =16; (bottom) 
standard deviation of PLIF intensity with transitional and turbulent 
locations marked. 


the single-shot images were examined. The flow was then 
labeled laminar, transitional, or turbulent based on the state of 
the flow at the most downstream location within the field of 
view of the camera. For laminar flows, xi am was recorded; for 
transitional flows, x trans was recorded. For turbulent flows, 
both x trans and x turb were recorded. The test matrix is plotted in 
Fig. 4, showing the range of Reynolds numbers and jet 
pressure ratios covered during the test. 

Most data were acquired with a 64 mm field of view. For 
selected test conditions, a wider field of view (165 mm) was 
employed in order to determine the jet behavior further 
downstream. A flow condition that appeared transitional 
within the smaller field of view might have become fully 
turbulent within the wider field of view. Test points are 
shown for both the 64 mm and 165 mm fields of view; 
therefore some test points are labeled with two labels (e.g. the 
test point for JPR=43 and unit Re exit = 0.6 x 10 6 m" 1 is labeled 
as both transitional and turbulent). 

In Fig. 4, points above a certain Reynolds number (about 3 
x 10 6 m" 1 ) appear to become fully turbulent, while points 
below a certain Reynolds number (about 0.3 x 10 6 m" 1 ) appear 
to remain laminar, which was the expected trend. However, 
this leaves a Reynolds number rang — spanning about an order 
of magnitude — over which flows exhibit mixed behavior. 
Clearly, additional parameters besides JPR and Re exit need to 
be considered, not only to predict if a flow will become 
turbulent, but also at what downstream distance the transition 
to turbulence will occur. 

These same data are plotted in Fig. 5 in terms of unit exit 
Reynolds number versus downstream distance, rather than jet 
pressure ratio. The general and expected trend appears — 
flows with lower unit Reynolds numbers are more likely to 
remain laminar. However, a clear relationship indicating the 
distance at which one might expect a given flow to undergo 
transition to turbulence is not readily apparent. For a many 
values of Re exit , transition occurred at different distances for 
different test conditions. For instance, some flows with a unit 
Reynolds numbers near 1 x 10 6 m -1 became fully turbulent in 
less than 20 mm, while others remained laminar beyond 160 


A A A 

♦ Laminar 

▲ A 

■ Transitional 

. 

A Turbulent 

▲ Aa ▲ 


■ ▲ 

A 

4 A 1 A 

A 

▲ m 


# 


A A ■ ■ 

♦ 

♦ 

♦ 

♦ 

tV 

A* 

a it 

♦ 


0 10 20 30 40 

Jet Pressure Ratio 


Figure 4. Unit exit Reynolds number versus jet pressure ratio. This plot 
shows the range of conditions over which these tests were conducted. A data 
point with a unit exit Reynolds number of nearly 15 and a JPR of 7.5 is not 

shown. 
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Figure 5. Unit exit Reynolds number versus distance downstream of the 
nozzle exit. 


mm. Note that laminar data fall at one of two distances; these 
correspond to the extent of the camera’s field of view, either 
64 mm or 165 mm. 

The lack of an obvious correlation in Fig. 5 led us to 
postulate that the jet behavior was a function not only of 
Reynolds number, but of jet pressure ratio as well. Scaling 
laws described by Yiiceil and Otiigen in Ref. [11] provide a 
means of redefining the Reynolds number to include the 
effects of JPR. Conditions at the nozzle exit serve as initial 
conditions for defining an “equivalent jet.” They assume that 
the initial jet is allowed to expand adiabatically until the 
pressure in the jet reaches the ambient pressure. At this 
point — which they call location 2 — all of the flow conditions 
are recalculated. Equation (7) can be combined with the 
equations for rescaled velocity V 2 , density p 2 , and temperature 
T 2 from Ref. [11] to arrive at an expression for the rescaled 
unit Reynolds number in terms of JPR and nozzle exit 
conditions. This expression is given in Eq. (9): 
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Fig. 6 shows the data plotted in terms of this rescaled unit 
Reynolds number, unit Re 2 . The relationship between this 
rescaled Reynolds number and transition length is more 
definitive. Independent power law fits to the transitional and 
turbulent data give R 2 values of better than 0.7. It is also 
important to note that all laminar points now lie below the 
fitted curve for turbulent test points. This indicates that, for a 
given Re 2 , there exists a distance beyond which the flow will 
likely be fully turbulent, and below which the flow will likely 
remain laminar. 

While attempting to find any systematic errors in the data, 
we reexamined the physical structures of the flows under 
investigation. We observed that flows with jet pressure ratios 
of less than about 3 exhibited different physical characteristics 


Figure 6. Correlation between rescaled unit Reynolds number and the 
distance downstream at which the jet flow became transitional (orange 
squares) and turbulent (red triangles). Power law fits to the data are shown. 
This plot includes data for jet pressure ratios ranging from 1.8 to 38. 



Figure 7. Single-shot PLIF images of underexpanded jet flows at five 
different jet pressure ratios. From top to bottom, these JPRs are: 36, 22, 8, 5, 
and 2. The presence of a Mach disk characterizes high JPR flows, while low 
JPR flows are characterized by a diamond shock pattern. Mach disks were 
observed in flows with JPRs greater than about 3. 



than flows with higher JPRs. Fig. 7 contains images of flows 
at five different jet pressure ratios. For the top four cases, 
flow structures including a barrel shock structure, a Mach disk 
and a shear layer can be seen. In the fourth image, the Mach 
disk is fairly small, and the flow is beginning to exhibit 
diamond shock features. In the bottom case, a Mach disk is no 
longer evident; the flow instead consists of a diamond shock 
pattern. 


When and if flows with high jet pressure ratios become 
turbulent, it appears that instabilities first develop in the shear 
layer, growing until they reach the core of the jet. One such 
case was illustrated in Fig. 3 for a flow with a JPR of 16. Fig. 
8 illustrates the effect that Reynolds number has on the 
location of the onset of transition for flows with JPRs of about 
16. The general trend is for flows at higher Reynolds numbers 
to undergo transition in a shorter distance. 



Figure 8. The effect of Reynolds number on transition to turbulence for a jet pressure ratio of about 16. Single-shot PLIF images (left) are displayed along with the 
corresponding standard deviation maps (right) for six values of unit Re exit (top to bottom): 0.62, 0.95, 1.8, 2.3, 2.7, and 14.8 x 10 6 m' 1 . 



Figure 9. PLIF image illustrating the “pinching off’ mode to turbulence, 
characteristic of flows distinguished by a diamond shock pattern. This flow 
has a JPR = 2. See Fig. 3 for an illustration of a turbulent flow at a JPR =16. 


Flows with low jet pressure ratios seemed to exhibit a 
different pattern when and if they became turbulent. Fig. 9 
illustrates this for a flow with a JPR of 2. As the flow 
becomes unstable, the diamond shock pattern appears to 
become distorted, and globules of jet fluid begin to be pinched 
off. 

Donaldson and Snedeker use a JPR of approximately 2 to 
delineate the boundary between moderately underexpanded 
and highly underexpanded jets [1]. They also note that the 
differences in physical flow structures — such as the 
appearance of a slip line and shear layer in highly 
underexpanded jets — affect jet diffusion, spreading, mixing, 
and dissipation rates. 

These physical differences between flows of different 
pressure ratios led us to consider only data points 
corresponding to JPRs of greater than 3 in our final 
correlation. The final correlation that was found between the 
rescaled unit Reynolds number and transition length is 
displayed in Fig. 10. The lower line shows the correlation 
between the rescaled Reynolds number and distance to 
transition; the upper line shows the distance to turbulence. 
Scatter in the data has been reduced by limiting the correlation 
to highly underexpanded jets with JPRs greater than 3. 

From these correlations, one can predict whether a given 
underexpanded jet will remain laminar over a certain distance, 
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Figure 10. Correlation between rescaled unit Reynolds number and the 
distance downstream at which the jet flow became transitional and turbulent. 
Data for jet pressure ratios less than 3 are not included here. 


or whether turbulent effects will become manifest. These 
predictions can provide valuable information to researchers 
running CFD simulations of such flows, indicating whether 
turbulence models need to be included in the simulation of a 
given flow. The correlations also help to define the range over 
which underexpanded jet flows will be neither fully laminar 
nor fully turbulent, regions where experimental measurements 
may prove to be most critical. 

IV. Conclusion 

In this work, we have demonstrated a powerful image 
processing technique to aid in the identification of turbulence 
using PLIF flow visualization images. This technique 
highlights turbulent regions of a flow by visualizing the 
standard deviation of fluorescence intensity throughout the 
flow. The image processing technique employed in this work 
could also readily be extended to studying transition and/or 
turbulence effects in other flows, particularly for tests in 
which spatially and temporally resolved flow visualization 
data is desired. In the current study, we have used the 
information gathered from such images to identify a 
correlation between a rescaled Reynolds number and 
downstream distance to transition for highly underexpanded 
axisymmetric sonic free jet flows. 
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